Mott transition and heavy fermion state in the pyrochlore Hubbard model 
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We investigate the interplay between geometrical frustration and strong electron correlation based 
upon the pyrochlore Hubbard model. In the half-filling case, using the perturbative expansion in 
terms of electron correlation, we show that the self-energy shows a divergent behavior leading the 
system into the Mott insulating state, in which quantum disordered spin liquid without magnetic 
long-range order realizes. In the hole-doped case, we obtain heavy-fermion-like Fermi liquid state. 
We also calculate the neutron cross section which is well consistent with recent neutron scattering 
experiments for itinerant pyrochlore systems. 
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I. INTRODUCTION 

Recently, the role played by geometrical frustration 
in strongly correlated electron systems has attracted 
renewed interests. EtQ The pyrochlore lattice, namely, 
corner-sharing tetrahedra (FIG.l), is one of the typi- 
cal system in which the geometrical frustration is crucial 
to determine its properties. Here we investigate effects 
of geometrical frustration on itinerant electron systems 
based upon the pyrochlore Hubbard model. One of the 
purpose of this paper is to investigate the Mott transi- 
tion in the absence of magnetic long-range order. Ac- 
cording to the study of the Heisenberg model on the py- 
rochlore lattice, it has been pointed out that the ground 
state of ±he insulating state is quantum disordered spin 
liquidllilj The critical character of the Mott transition 
without magnetic order is quite different from that ac- 
companying antiferromagnetic one. Moreover the char- 
acter of the quantum spin liquid has not yet been unveiled 
sufficiently It is expected that the study from electron 
systems will shed new light on this subject. 

Another purpose of this paper is to reveal how the ge- 
ometrical frustration affects strong electron correlation 
effects in the metallic state. Some recent experiments on 
the pyrochlore itinerant electron systems have reported 
that remarkable heavy-fermion-like behaviors manifest in 
these systems, implying the crucial role of the geometrical 
frustration.Ell3 Here we show that the heavy-fermion state 
is realized in the pyrochlore Hubbard model in the vicin- 
ity of the half-filling. We also show that the dynamical 



spin susceptibility obtained from our model away from 
the half-filling is consistent with recent neutron scatter- 
ing experiments for itinerant pyrochlore systems .c3 

The organization of this paper is as follows. The model 
and the basic method are given in Sec. II. In Sec. Ill, we 
discuss about the Mott transition and the realization of 
the quantum disordered spin liquid in the half-filling case. 
In Sec. IV, we consider the hole-doped case focusing on 
the heavy fermion state and its magnetic properties. 




FIG. 1. Pyrochlore lattice. 



II. MODEL HAMILTONIAN AND METHOD 



The Hamiltonian of the model is given by, 
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x yJ cos(k y ~k z ) cos(k x ~ k z ) 

„r> ,1 I cos(k x + k y ) cos(k x + k z ) cos(k y + k z ) 

Vka I CO s(fcj, — k z ) cos(k x + k z ) cos{k x — ky) 

cos(k x — k z ) cos(k y + k z ) cos(k x — k y ) 



(1) 



Here Ci UrT {<y iva ) is an annihilation (creation) operator 
of electrons with spin a on the f-th site of the i-th tetra- 
hedron. = ( c m ct > c L CT ' c !,3<t> c L4 CT )- Thc kinctic 

term is diagonalized as, H kin = ^ fco - Y,t=i E kua\ wa a kv< j, 



w here E kl = E k2 = 1, E k3 = -1 + y/l + tk, E k4 = -1 - 
y/1 + t k with t k — cos(2k x ) cos(2k y ) +cos(2k y ) cos(2k z ) + 
cos(2fc z ) cos(2fca;). The basis in the diagonalized space 
is obtained from the canonical transformation, ipt = 
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*51v=iSv{k)a kv(J . Using the abbreviation, s(x ± y) = 
sm(k x ± k y ) etc., we write down s\ 2 (k), 

s 1 (k) = (s(x + z), s(y - z), -s(x + y),0)/n + , (2) 

where n± = y/ s(x ± z) 2 + s(y =F z) 2 + s(x + y) 2 , and, 

s 2 (k) = (s(x + z)s(x — z)s(y — z) 

-s(y + z)(s(y - zf + s(x + y) 2 ), 

s(x + z)s(y + z)s(y — z) 

—six — z)(s(x + z) 2 + s(x + y) 2 ), 

— s(x + y)(s(x + z)s(y + z) + s(x — z)s(y — z)), 

s(x + y)n 2 + )/n 2 , (3) 

for k x + k y 7^ 0, where n 2 = n + [n' 2 yn 2 _ — (s(x + z)s(y + 
z) + s(x — z)s(y — z)) 2 ]^ , and, 

s 2 (k) = (-s(2x),-s(2x),2s(x - z),2s(x + z)) 

x l/ v / 2s(2x) 2 + 4s(x - z) 2 + As{x + z) 2 (4) 

for k x + k y = 0. The expressions of s%{k) and Si{k) are 
very complicated. However in the following, we need only 
s^(k) for small fc, which is given by, 

s 3(k) = (—k x — ky + k z , k x + ky + k z , 

kx ~l~ ky k x ky fc 2 )/2|fc|. (5) 

At the half-filling, in the absence of electron-electron 
interaction, Ek3 and Eki are filled completely, and the 
two degenerate flat bands are empty. Ek3 touches with 
the two flat bands at the T point of the Brillouin zone. 
This state is a gapless band insulator. We now consider 
the effect of electron interaction on this state. Since, in 
the case of the half-filling or in the vicinity of the half- 
filling, the band E^a is sufficiently far from the Fermi 
level, it does not affect low energy properties and is negli- 
gible in the following argument. Then, the single-particle 
Green's functions of a^o-electron are determined by the 
following equation. 

+ n- E kX )S, x - Ej> x >( £ )]Gf > (e) - V (6) 

A=l 

We put the chemical potential \i = 1. In general, the off- 
diagonal self-energy E^^/i ^ v) is not negligible. How- 
ever, in the vicinity of the T point, at which the most 
important scattering processes occur, the off-diagonal 
terms vanish because of the momentum dependence of 
s v {k). Thus in the following argument, we neglect the 
off-diagonal self-energy. 

III. THE HALF-FILLING CASE 

A. Self-energy 

To investigate the electronic state at the half-filling, 
we calculate the self-energy. Recently, Isoda and Mori 



obtained the diagonal self-energy of this model up to the 
second order in C/J13 which is given by, ReE(e) ~ \J~-e 
for e < 0. However, as will be shown here, the higher 
order corrections to the self-energy give rise more singu- 
lar contributions and a drastic change of the electronic 
state, namely, the Mott transition. 
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FIG. 2. The most singular diagrams of the self-energy 
up to the forth-order in U. a denotes a3-electron. b 
denotes a\ or a2-electron. The dotted line denotes the 
interaction U. 

At the half-filling, the Ek3-ba,nd is completely filled, 
while the flat bands are empty. Hence, at zero temper- 
ature, particle-particle channel between a3-electron and 
ai or a2-electron vanishes. On the other hand, particle- 
hole channel gives non-vanishing contributions only for 
the pair of a3-electron and a\ or Q,2-electron. Then the 
most singular diagrams up to the forth order in U are 
those shown in FIG. 2. Since the Fermi level is at the 
position between the two flat bands and i?fc3-band, we 
shift the chemical potential infinitesimally, fj, — * 1 — S. 
After integrating over energy and momentum, we take 
the limit 5 — * 0. This procedure is required in this singu- 
lar perturbative calculation to fix the density of electrons 
n = 1. Neglecting the momentum dependence of s v (k) 
does not affect the leading singular behavipis. Then, we 
can carry out the calculation analytically.Ea For e < 0, 
we obtain, 

. JJ3 Jji 

ReE(e) ~ const. - cU 2 ^\ ~ 0—^= + c— - T , (7) 

where c is a positive constant of order unity. The real 
part of the self-energy for hole-like excitations is diver- 
gent. The same divergent behavior appears in the imag- 
inary part, however, for particle-like excitations. Higher 
order corrections also show stronger divergence. This 
divergent behavior means that the unperturbed state is 
unstable in the presence of electron correlation and that 
the single particle energy gap is generated. To see this we 
take only up to the third order term of the self-energy. 
Hole-like excitations are possible only for i^-band in 
the half-filling case at zero temperature. Thus the sin- 
gle particle energy of this band measured from the Fermi 
level is changed to ~ — U 2 = —A for k = 0. 
band is pushed down to the lower energy, and the energy 
gap opens between the flat bands and -Efc3-band. This 
gap generation due to electron correlation signifies that 
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the systems is in the Mott insulating state. Although the 
exact magnitude of the gap should be determined by tak- 
ing into account higher order corrections in U, the gapful 
state is the self-consistent solution for non-vanishing U. 
Thus an infinitesimally small U drives the system into the 
Mott insulator at the half-filling. In the ground state, the 
two bands below the gap, Ek3, -Em, are filled with up and 
down spins. Hence the system is in the spin singlet state. 
It should be stressed that, after the mass gap generation, 
the singularities of the perturbative expansion are elimi- 
nated, and as a result, the gapful state is stabilized even 
if one take into account higher order corrections, which 
just renormalize the magnitude of the gap. 



B. Quantum spin liquid state 



As the electron filling approaches the half-filling value 
\fi\ — > 0, the mass enhancement factor shows the diver- 
gent behavior indicating the precursor of the Mott transi- 
tion. Such a large mass enhancement is actually observed 
in the specific heat measurement for some pyrochlore itin- 
erant systems like Y(Sc)Mn2.El 




FIG. 3. The most singular diagrams of spin-spin corre- 
lation functions up to the third-order in U. 



To see the magnetic property of the Mott in- 
sulating state obtained in the previous section, we 
compute spin-spin correlation functions, \^(q,Lu) — 

(S^(q,u)S-(-q,-u)}, where S+ = J2 k c l+ q ^ c ^l- We 
carry out perturbative expansion in terms of U using 
the Green's function of the gapful state, G^ 33 ^(e) ~ 
[s — Ek3 + A] -1 , as the unperturbed propagator. Calcu- 
lating the diagrams up to the third order in U, we have, 
Imx(q, k>) ~ y/u — AO(w-A) with Q(x) a step function. 
The spin gap exists between the spin singlet ground state 
and the triplet state. Thus the system is in the quan- 
tum disordered spin liquid state without magnetic long- 
range order. The result may not be changed qualitatively 
by higher order corrections. These results are consis- 
tent with the-prpdous studies for the pyrochlore Heisen- 
berg model.aE3il3 The singularity of the self-energy (0) 
is deeply related with the accidental degeneracy at the T 
point in the momentum space, which may cause an insta- 
bility of lattice structure. This property may be relevant 
to the metal- insulator transitioipaf TI2RU2O7, which ac- 
companies a structure change.BEa 



IV. THE HOLE-DOLED CASE 

A. Heavy fermion state 

In the hole-doped case, the singularities appeared in 
the half-filling case are eliminated by the presence of the 
cutoff which is the chemical potential \i measured from 
the half-filling level. Then the Fermi liquid metallic state 
is realized. However, in the vicinity of the half-filling, the 
effective mass is much enhanced by electron correlation. 
Up to the third order in U, the leading term of the mass 
enhancement factor is, 

dReZje) + (g) 
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FIG. 4. (a) The structure factor 312(9) plotted in the 
[001]- [110] plane. The black line denotes the zone bound- 
ary. The brighter regions have the stronger intensity, (b) 
The neutron cross section. 
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B. Magnetic properties and the neutron cross 
section 

The largely enhanced effective mass obtained in the 
previous section is a remarkable property of geometri- 
cally frustrated electron systems in which several modes 
of spin fluctuations compete with each other, and mag- 
netic long-range order is suppressed. The spin fluctua- 
tion is almost localized in real space as in the case of 
/-electron based heavy fermion systems. This property 
is seen more clearly in the momentum dependence of the 
spin-spin correlation function. We calculate it by the per- 
turbative calculation. Since the magnetic frustration is 
so strong, any approximations which neglect the coupling 
between several modes of spin fluctuation like random 
phase approximation(RPA) or fluctuation exchange ap- 
proximation(FLEX) are invalid for this system. Actually 
the most singular contributions up to the order of 0(U 3 ) 
comes from the diagrams shown in FIG. 3. To simplify the 
calculation we exploit the following approximation. Near 
the half-filling, the momenta running in G^^-line take 
small values in the vicinity of the Fermi surface. More- 
over for small k, S3 (A;) does not depend on |fc|. Thus we 
can separate the integral over |fc| and the angle for the 
momenta running in G( 33 '-line. Then, the spin-spin cor- 
relation function is decomposed into the momentum part 
and the energy part, Xnv(Q) u ) = fi^dn^io)- Such a fac- 
torization of the spin-spin correlation function is indeed 
observed in an_inelastic neutron scattering experiment 
for Y(Sc)Mn2.li3 The energy part is given by, 



/(w)~-aitf 2 log 



2/z| 



E c 



-a 2 U 3 



1 



LO - 



2[i 



(9) 



for small u>. Here 01, a 2 are constants. Thus in the vicin- 
ity of the half-filling, spin fluctuation is much enhanced. 
The presence of the giant spin fluctuation is also con- 
sistent with, several experiments for itinerant pyrochlore 
systems I!rlii3 It is noted that such a large enhancement 
of X/j,v(<1,v) is not obtained by RPA or FLEX approx- 
imations for the appropriate value of U. We calculate 
the structure factor g^vio) numerically. In FIG. 4(a), we 
plot the intensity of 512(9) in the [110]-[001] plane of the 
reciprocal lattice. 512(9) is almost constant in the [001]- 
direction, indicating that the spin fluctuation is strongly 
localized in this direction. It also has a small peak on the 
9.t = Qy = line showing the presence of small fluctuation 
toward a collinear magnetic order, as was pointed out in 
the study of the Heisenberg model.cl However this mode 
competes with the other modes which exist in 313(g), 
523(9), 514(9), 524(9), and thus such a magnetic order is 
suppressed. 

In FIG. 4(b), we show the calculated result of the neu- 
tron cross section, £^ v e~ lGtyR ^~ R ^x^v{l^ w ), m the 
[110]- [001] plane. Here G is a reciprocal lattice vec- 
tor. The positions of the maximum are almost consistent 
with the neutron scattering experiment for Y(Sc)Mn2.t£l 
Note that this structure is very similar to that found for 



the Heisenberg model.cl Thus the itineracy of electron af- 
fect little the structure in the momentum space, though 
the dynamical properties between the Mott insulating 
state and the metallic state are quite different; namely, 
the former is gapful and the latter is gapless. The mo- 
mentum dependence of X^u(q,oj) implies that the spin 
fluctuation is localized in a tetrahedron forming collec- 
tive 4-spin singlet .Ej'S Such a localized character of spin 
fluctuation is very similar to magnetic properties of /- 
electron based heavy fermion systems. However, in the 
latter systems, magnetic long-range order is suppressed 
by the presence of Kondo temperature higher than mag- 
netic ordering temperature. In our system, the geomet- 
rical frustration is crucial. 



C. Compressibility 



Here we discuss the charge response of the heavy- 
fcrmion state near the half-filling. In general, the charge 
susceptibility Xc in the vicinity of the Mott metal- 
insulator transition point shows quite different behav- 
iors depending on the density of states near the Fermi 
levelElJ In our systems, the bare density of states in the 
hole-doped case is non-singular, though four-point ver- 
tices are much enhanced by electron correlation. The 
perturbative calculation of the charge susceptibility gives 
singular terms in the limit of /1 — * 0, namely, Xc ~ 
d U 2 log I fi I + C 2 U 3 1 1 n I + • • • . Thus it is highly non- trivial 
how the charge susceptibility behaves, as the electron fill- 
ing approaches the half-filling value. To see this we resum 
the most singular terms of the perturbative expansion in 
U. As was done before, we neglect the momentum depen- 
dence of s„. We first consider the terms which contain 
only one irreducible four-point vertex. The most singu- 
lar terms of this type are diagrammatically expressed as 
FIG. 5 (a) These terms are calculated as, 



= c 2 u 2 log + c 3 [/ 3 /M + c 4 t/7M 2 



(10) 



Here C2,C3, and C4 are positive constants. On the 
other hand, the term which contains n irreducible four- 
point vertices shown in FIG. 5(b) gives the contribution 

~ [ — cXcTVImI^ - with c a positive constant. Summing 
up all most singular terms, we have, 



X" 



(11) 



Thus as the electron filling approaches the half-filling 
value /j, — ► 0, the charge susceptibility decreases toward 
zero, Xc —> 0, indicating that the system becomes incom- 
pressible. 
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FIG. 5. Diagrams of the charge susceptibility. Ti rr is 
the irreducible 4-point vertex. 



V. SUMMARY 

In this paper, we investigate the Mott transition, 
strong correlation effects and spin dynamics of the py- 
rochlore Hubbard model. We show that, at the half- 
filling, electron-electron interaction leads the system into 
the Mott insulator in which the quantum disordered spin 
liquid state realizes. In the hole-doped case, we show that 
the effective mass is anomalously enhanced in the vicin- 
ity of the half-filling signifying the heavy fermion state. 
We have also calculated the neutron cross sectin of this 
state which is consistent with the recent experiment for 
some itinerant pyrochlore systems like Y(Sc)Mn2. 

Although we carried out the perturbative calculation 
up to the fourth order in U , the singular divergent behav- 
ior of the self-energy shown in Sec. Ill implies that we 
need to resum all order singular diagrams. Since the most 
singular diagrams have the specific structure as discussed 
in Sec. Ill, we can carry out this program in principle. 
We would like to address this issue in the near future. 
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